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Precise determination of the strangeness magnetic moment of the nucleon 
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By combining the constraints of charge symmetry with new chiral extrapolation techniques and 
recent low mass lattice QCD simulations of the individual quark contributions to the magnetic 
moments of the nucleon octet, we obtain a precise determination of the strange magnetic moment of 
the proton. The result, namely G\.i — —0.046 ±0.019 /i]v, is consistent with the latest experimental 
measurements but an order of magnitude more precise. This poses a tremendous challenge for future 
experiments. 



There is currently enormous interest in the determina- 
tion of the strangeness content of the nucleon. It is cru- 
cial to our understanding of QCD to determine precisely 
the role played by heavier, non-valence flavors. On the 
experimental side new results on strangeness in the nu- 
cleon have been reported recently from JLab (HAPPEX) 
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cieon nave oeen reportea recently irom Ji^aD (^MArriijAj 
[3 and MIT-Bates (SAMPLE) [2. In the near future we 
can expect even more precise results from the A4 experi- 
ment at Mainz as well as GO and HAPPEX2 at JLab. By 
contrast, the theoretical situation is somewhat confused, 
with the predictions of various quark models covering an 
enormous range. Direct calculations within lattice QCD 
have not yet helped to clarify the situation, with values 
for G%i ranging from -0.28 ± 0.10 to -h0.05 ± 0.06 0. 

We take a different approach, exploiting the ad- 
vances in lattice QCD which have enabled quenched 
QCD (QQCD) simulations of magnetic moments at pion 
masses as low as 0.3-0.4 GeV IMHBIgl) as well as the de- 
velopment of new chiral extrapolation techniques J9|, Il0l| . 
Using these techniques we determine, in full QCD, the 
ratios of the valence u-quark contribution to the mag- 
netic moment of the physical proton to that in the S'^ 
and of the valence u quark in the physical neutron to that 
in the 'EP . From these ratios, the experimental values of 
the octet moments and charge symmetry we deduce a 
new theoretical value for G%j which is precise - setting 
a tremendous challenge for the next generation of parity 
violation experiments. 

As illustrated in Fig. Q] the three point function re- 
quired to extract a magnetic moment in lattice QCD in- 
volves two topologically distinct processes. (Of course, in 
full QCD these diagrams incorporate an arbitrary num- 
ber of gluons and quark loops.) The left-hand diagram 
illustrates the connected insertion of the current to one 
of the "valence" quarks of the baryon. In the right-hand 
diagram the external field couples to a quark loop. The 
latter process, where the loop involves an s-quark, is en- 
tirely responsible for G|j. 

Under the assumption of charge symmetry |lll |. the 



FIG. 1: Diagrams illustrating the two topologically different 
insertions of the current within the framework of lattice QCD. 

magnetic moments of the octet baryons satisfy [lj| : 
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Here, p and S^ arc the physical magnetic moments of 
the proton and S^, and similarly for the other baryons. 
The valence u-quark sector magnetic moment in the pro- 
ton, corresponding to the LHS of Fig. ^ is denoted 
uP. Charge symmetry has been used to replace the d- 
quark contribution in the neutron by u^, d in the E~ 
by u in the S^ ( w^), and so on. The labels on quark 
magnetic moments allow for the environment sensitiv- 
ity implicit in the three-point function 12. 13] . That is, 
the naive expectations of the constituent quark model, 
namely u^ /vP = uP" ju" = 1, may not be satisfied. The 
total contribution from quark- loops. On, contains sea- 
quark-loop contributions (RHS of Fig. ^ from u, d and 
s quarks. By definition 
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where the ratio of s- and d-quark loops, ^R^ = 
^G^j/^G^j, is expected to lie in the range (0,1). Note 
that, in deriving Eq.Q, we have used charge symmetry 



to set ^G'^j = ^Gj^j. Since the chiral coefRcients for the d 
and s loops in the RHS of Fig. 1 are identical, the main 
difference comes from the mass of the K compared with 
that of the tt. 

With a little algebra On, and hence G\^{= ^'G\j), may 
be isolated from Eqs. ^ and ||2J): 
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Incorporating the experimentally measured baryon mo- 
ments |l4| (in nuclear magnetons, /iAr), Eqs. Q and (Q 
become: 
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sequences of QCD, under the assumption of charge sym- 
metry. 

Equating JB)) and ^ provides a linear relationship be- 
tween uP /vP and vP" ju" which must be satisfied within 
QCD under the assumption of charge symmetry. Figure 
12 displays this relationship by the dashed and solid line. 
Since this line does not pass through the point (1.0, 1.0), 
corresponding to the simple quark model assumption of 
universality, there must be an environment effect exceed- 
ing 12% in both ratios or approaching 20% or more in at 
least one of the ratios. Indeed, a positive value for 01^(0) 
would require an environment sensitivity exceeding 70% 
in the w"/u" ratio! 

The numerical simulations of the electromagnetic form 
factors presented here are carried out using the Fat Link 
Irrelevant Clover (FLIC) fermion action 'j?, 'fq in which 
the irrelevant operators, introduced to remove fermion 
doublers and lattice spacing artifacts, are constructed 
with APE smeared links (151] . Perturbative renormal- 
izations are small for smeared links and the mean-field 
improved coefficients used here are sufficient to remove 
0(a) errors from the lattice fermion action |l6| . 

The 0(a)-improved conserved vector current [13 is 
used. Nonperturbative improvement is achieved via the 
FLIC procedure, where the terms of the Noether cur- 
rent having their origin in the irrelevant operators of the 
fermion action are constructed with mean-field improved 
APE smeared links. The results presented here are ob- 
tained using established techniques [l3| from a sample of 
400 20^ X 40 mean-field 0(a^)-improved Luscher-Weisz 
[igj l gauge field configurations having a lattice spacing of 
0.128 fm, determined by the Sommer scale ro = 0.50 fm. 

One of the major challenges in connecting lattice cal- 
culations of hadronic properties with the physical world 
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FIG. 2: The constraint (dashed Gm(0) < 0, solid Gm(0) > 0) 
on the ratios vF juP and u^ /u~ implied by charge symmetry 
and experimental moments. Experimental uncertainties are 
indicated by the dotted bounds. The assumption of environ- 
ment independent quark moments is indicated by the crossed 
square. Our final result (chiral corrected extrapolation of lat- 
tice results) is illustrated by the filled square on the charge 
symmetry line. 

[20| is that currently accessible quark masses are much 
larger than the physical values. Our present analysis has 
been made possible by a significant breakthrough in the 
regularization of the chiral loop contributions to hadron 
observables |9|, llOl |21| . Through the process of regulat- 
ing loop integrals via a finite-range regulator 9, 22], the 
chiral expansion is effectively re-summed to produce an 
expansion with vastly improved convergence properties. 
In particular, we extrapolate FLIC fermion calculations 
of the valence quark contributions to baryon moments 
(m^, m", m^, u") to the physical mass regime. We select 
the dipole-vertex FRR with A = 0.8 GeV, which yields 
the best simultaneous description of both quenched and 
dynamical simulation results J23l | . 

Separation of the valence and sea-quark-loop contri- 
butions to the meson cloud of full QCD hadrons is a 
non-trivial task. We use the diagrammatic method for 
evaluating the quenched chiral coefficients of leading non- 
analytic terms in heavy-baryon quenched xPT (Q^PT ) 
|2J, I23 • The valence contributions (key to this analysis) 
are obtained by removing the direct-current coupling to 
sea-quark loops from the total contributions. Upon fur- 
ther removing "indirect sea-quark loop" contributions, 
where a valence quark forms a meson composed with a 
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FIG. 3: Diagrams providing the leading contributions to 
the chiral expansion of baryon magnetic moments (upper di- 
agrams) and their associated quark flows (lower diagrams) in 
QQCD. 
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FIG. 4: The contribution of a single u quark (with unit 
charge) to the magnetic moment of the proton. Lattice sim- 
ulation results (square symbols for rn^ > 0.05 GeV) are ex- 
trapolated to the physical point (vertical dashed line) in finite- 
volume QQCD as well as infinite volume QQCD, valence and 
full QCD - see text for details. Extrapolated values at the 
physical pion mass (vertical dashed line), are offset for clarity. 




FIG. 5: Correcting QxPT (upper) to the valence sector of 
full QCD (lower diagrams). We remove quenched negative- 
metric rj' contributions and adjust the chiral coefficients of vr 
and K loops to account for the coupling of a valence quark to 
the photon in a meson made from a sea-quark loop. Coupling 
to the anti-quark in the bottom-right diagram is also included 
in the valence contribution of full QCD. 

sea-quark loop, one obtains the "quenched valence" con- 
tributions - the conventional view of the quenched ap- 
proximation. 

Figure O displays the diagrams providing the leading 
contributions to the chiral expansion of baryon magnetic 
moments (upper diagrams) and their associated quark 
flows in quenched QCD (QQCD). The associated chiral 
expansion for the proton magnetic moment, ^p, has the 
form 

Mp = flg +^pX»j'-f.,'(nT-7r, A) +X7rB-fs(m-7r, A) 

+XKBlB{mK,A) +a2ml + a^ml. (8) 

where the repeated index, B, sums over allowed baryon 
octet and decuplet intermediate states. Loop integrals 
denoted by / are defined by 
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FIG. 6: The contribution of a single u quark (with unit 
charge) to the magnetic moment of the E^ . Curves and sym- 
bols are as for Fig. 2] 

where Abn is the relevant baryon mass-splitting and 
u{k,A) is the dipole-vertex regulator. The coefficients, 
X, denote the known model-independent coefficients of 
the LNA term for tt and K mesons 25, 261. We take 
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and use the physical values to define 

The m^ term in Eq. (jHl allows for some curvature 
associated with the Dirac moment of the baryon, which 
should go as l/m^ for moderately large quark masses. 

Figure H illustrates a fit of FRR QxPT to the FLIC 
fermion lattice results (solid curve), where only the dis- 
crete momenta allowed in the finite volume of the lattice 
are summed in performing the loop integral. The long- 
dashed curve that also runs through the lattice results 
corresponds to replacing the discrete momentum sum by 
the infinite-volume, continuous momentum integral. For 
all but the lightest quark mass, finite volume effects are 
negligible. 

The coefficients of the residual expansion, Oq , a2, a^, 
show excellent signs of convergence. For example, the fit 
to u^ yields values 1.48(7), -0.90(23), and 0.42(19) in 
appropriate powers of GeV, respectively. Incorporating 
baryon mass splittings into the kaon loop contributions 
is essential - e.g., the contribution of S — > NK is almost 
doubled when the "E — N mass splitting is included. 

FigureElillustrates the considerations in correcting the 
quenched u-quark contribution to yield the valence u- 
quark contribution in full QCD. The removal of quenched 
rj' contributions and the appropriate adjustment of n and 
K loop coefficients |2J, |2^ |2a| provides the dot-dash 
curve of Fig. 0] This is our best estimate of the valence 
M-quark contribution (connected insertion) to the proton 
magnetic moment of full QCD. Finally, the disconnected 
insertion of the current is included to estimate the total 
contribution of the M-quark sector to the proton mag- 
netic moment |2J,|23,|2y| (fine dash-dot curve in Fig.0J. 
Figure El displays similar results for the E+ . 

From these chiral extrapolations, we estimate the ra- 
tios of the valence (connected) w-quark contributions. 
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FIG. 7: The one standard deviation agreement between 
the FRRxPT corrected lattice simulation results (■) and 
the experimentally measured baryon magnetic moments (•). 
Quenched (♦) and finite- volume quenched (H) results are also 
illustrated. 

uP ji^ and u^ jvr . The final results 

^T = 1.092 ±0.030 and ^ = 1.254 ± 0.124 (11) 

are plotted in Fig. [3 The precision of these results follows 
from the use of correlated ratios of moments which act to 
reduce uncertainties associated with the lattice spacing, 
the regulator mass and statistical fluctuations [23|. This 
result leaves no doubt that G^ is negative. The fact that 
this point lies exactly on the constraint curve is highly 
nontrivial, and provides a robust check of the validity of 
the analysis techniques presented here. 

As a further check, in Fig. [7| we compare the lattice 
QCD predictions of the baryon magnetic moments con- 
structed from chirally-corrected extrapolations of the in- 
dividual quark sectors. The results display an unprece- 
dented level of agreement with experiment. We note that 
the experimental constraints on liP and u" emphasized 
by Wong !23| are both satisfied precisely. 

While G%i is most certainly negative, it remains to 
set the magnitude. This requires an estimate of the 
strange to light sea-quark loop contributions, ^R% Ear- 
lier estimates of ^i?^ were based on the constituent quark 
model. A more reliable approach is to estimate the loops 
using the same successful model invoked to correct the 
quenched results to full QCD [liJ, |23| as illustrated in 
Fig. H Allowing the dipole mass parameter to vary 
between 0.6 and 1.0 GeV provides ^i?^ = Glj/Gjj = 
0.139 ± 0.042. A complete analysis of the errors asso- 
ciated with the determination of G\,j using Eqs. I@J, (0) 
and (|ll(l is reported in Ref. ;27] . The uncertainty is dom- 
inated by the statistical errors included in Eq. (|ll|l and 
the uncertainty just noted for ^i?^. The final result for 
the strangeness magnetic moment of the nucleon is 



This precise value sets a tremendous challenge for the 
next generation of parity violation experiments. 



We thank the Australian Partnership for Advanced 
Computing (APAC) for generous grants of supercom- 
puter time which have enabled this project. This work 
is supported by the Australian Research Council and by 
DOE contract DE-AC05-84ER40150, under which SURA 
operates Jefferson Laboratory. 



Gli = -0.046 ± 0.019 ^1N 
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